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The Applicability of PQCD and NLO Power Corrections for Pion Form Factor 

Tsung-Wen Yeh* 
Institute of Physics, National Chiao- Tung University, Hsinchu, Taiwan 300, Republic of China 

In many years ago, Isgur and Llewellyn Smith addressed that PQCD is inapplicable to exclusive 
processes pJ,01, such as the pion form factor. The inain problem is that the asymptotic of PQCD 



is only about one fourth of the experimental value |3l-n]. We reexamine this PQCD deep problem. 
By including NLO power corrections to the pion form factor, we may arrive at a perturbative ex- 
planation for the data. The key realization is that we need to interprete that the strong interaction 
coupling constant involved in the PQCD result should be taken as an effective coupling constant 
under nonperturbative QCD vaccum. This implies that one can equally identify the relevant scale 
for the effective coupling constant as the factorization scale about 1 GeV. We also find that the 
average momentum fraction variable locates about 0.5, which is in favor of the asymptotic pion 
wave function. By employing photon-pion form factor with NLO power corrections to factorize out 
the nonperturbative effects involved in the effective coupling constant, we can extract an effective 
?C ' running coupling constant, which represents an effective coupling involving in the hard scattering 

^vj , subprocesses. The difference between the effective running coupling constant and the usual pertur- 

bative running coupling constant (Aqqj^ = 0.3 GeV) is very small for Q^ > I GeV^. The effective 
^ , running coupling constant Os/n is smaller than 0.2 for Q^ > I GeV'^. This directly showes that 

PQCD is applicable to exclusive processes at energy Q^ > 1 GeV^. In summary, with NLO power 
0^ I corrections, PQCD can completely explain the Q^ spectrum of pion form factor. 
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I. INTRODUCTION 



ri . The exclusive process plays an important role in understanding the strong interactions. The detail analysis of 

p) i' exclusive processes not only can reveal the constituents of hadron but also can untangle the underlying dynamics. 
' , One important progress in this respect is the prediction of perturbative QCD (PQCD) for exclusive process in high 
br^' energy limit MM- In PQCD formula, the transition amplitude can be expressed as a convolution integral over initial 
r^ , and final state meson wave fucntions and a hard funtion. The hard function can be perturbatively calculable, while 
the meson wave functions can only be determined, nonperturbatively. 

The pion form factor has been investigated in the framework of PQCD |q-[l7||. In the experimentally accessable 
energy region of few GeV^, the asymptotic of PQCD is only about one fourth of the experimental value. This leads 
^ ' to two problems for PQCD. One problem relates to the self-consistency of PQCD and the other problem requires 
increasing the maginutde of PQCD. The former was first addressed by Isgur and Llywellyn Smith |Q,p| that the most 
contributions come from the soft end-point regions or small values of the fraction variables, indicating the perturbative 
calculation inconsistent. This problem was partially resolved by Li and Sterman S who proposed a modified PQCD 
formulation, in which the physics of the end-point regions are carefully dealt with and the associated contributions 
are suppressed by a resummation over soft radiative corrections, the Sudakov factor. It was found that, with the 
transverse degrees of freedom playing the role of infrared cut-off, the PQCD contribution becomes self-consistent for 
momentum transfers as low as few GeV. However, note that the magnitude of modified PQCD is still smaller than 
the experimental value by a factor of 2 ~ 4. 

The general wisdom to increase the magnitude of PQCD is by invoking the power correction [J. However, the 
conclusion of past results was rather pessimistic. Even including twist-3 power corrections, it is still difficult to 
accomodate the data. More worsely, introducing the twist-3 contribtion imposes another problem that the twist-3 
contributions may exceed the leading twist contribution at intermediate energy region of Q^ being 2 ~ 15 GeV^ 
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'^The contributions come from the NLO in as is 20% ~ 30% 



|g2|-Q or 2 ^ 40 GeV^ (for the parton carrying transverse momenta) |T^. This is because there involves a huge 
factor rn^/mo « 1000 for the twist-3 contribution, where m,r is the pion mass and ttiq is the average constituent quark 
masss of the pion. The strange behavior of the twist-3 contribution violates the principle of perturbation theory: the 
magnitude of subleading terms should be smaller than that of the leading ones. 

The reason why the twist-3 contribution is inconsistent with perturbation theory may reside in the method of calcu- 
lation. In the previous approaches [|l2|- |l4| , |l7| , the spin projections 75 and 7517^'' are employed and the normalizations 
of related distribution amplitudes (DAs) are determined with the help of chiral perturbation theory. However, these 
methods gave no garantees on factorization of hard and soft radiative corrections for a twist-3 contribution. The 
failure of these methods in explaining the data may indicate the incorrectness of factorization. 

In a recent power expansion method developed by us [|l8| , it showes that the calculation of high twist contribution 
is compatible with QCD factorization. There requires more perturbative contributions than spin projections. The 
most important feature of this method is that the partonic interpretation for high twist contribtion is retained. For 
the pion form factor, the NLO power corrections calculated by this method have reverse signs and inverse magitudes 
against to those results of the previous mentioned methods. These differences give us a opportunity to explain the 
data by only employing PQCD. It is noted that the NLO contribution derived in this text is smaller than the leading 
order contribution over whole energy range of the experiment. 

Our organization is as follows. We investigate power expansion and QCD factorization for process 7*7r — > tt in 
Sec. II. This is a preliminary for Sec. III. The pion form factor Ft^{Q'^) upto order of 0{Q~'^) are evaluated in Sec. III. 
The proof showing that PQCD is applicable to pion form factor is present in here. A discussion for the reason why 
Isgur and Llewellyn Smith made their conclusion is explained in Sec. IV. Sec. V is devoted to conclusions. 

II. COLLINEAR EXPANSION AND QCD FACTORIZATION 

In this section, we will describe our approach of power expansion for 7*7r — > tt. The method we employ is called the 
collinear expansion |18|-p0[. In the following, we first sketch the procedures of how to pcrfom the coUinear expansion 
for 7*7r —> TT. Then, we will show the collinear expansion is compatible with the QCD factorization, which is important 
for giving a reliable perturbative calculation in QCD. 

A. Tree Level Collinear Expansion 

Let a = (f)* {P2, k2) ^ ap{ki, k2) ® (j>{Pi, ki) represent the lowest order amplitude for j*{q)Tr{Pi) — > tt{P2) as depicted 
in Fig. 1(a). The (Tp(fci, fc2) denotes the amplitude for partonic subprocess and the <t>{ki), i = 1,2, represents the pion 
DAs. The (S means convolution integral over the loop momenta fc^. To pick out the leading contribution, we assign 
the momenta of the initial and final state pions in the following way. We let Pi = Qvi where the dimensionless vectors 
obey vf = 0, vi ■ V2 = ^{l — cos 9). 9 is the center of mass scattering angle. The internal loop momenta k'^, i — 1,2 
are parameterized as 

k^ = XiQv^ + yiv[ + kij_ (1) 

Xi are dimensionless number of order unit, while yi have dimensions of mass. The vector u- are in the direction of the 
opposite- moving external vectors such that Vi ■ v'^ = 1, {vl)'^ = and v[ ■ vi^ = vi ■ V2, and also eijVi • u' = ^(1-1- cos 9) 
for i,j — 1,2. yi are solved by keeping kf invariant 

y,^ ^ ^i + ^'-L . (2) 

ZiXi^o^ 

The first step is to perform Taylor expansion for the parton amplitude 

<Jp{ki) = 'Spih = XiPi) + (CTp)a(xi, Xi)wf^,kf H (3) 

where we have assumed the low energy theorem 



dk 



a ^'^pK'^i 



= i(Jp)aiXi,Xi) . (4) 

ki=XiPi 



and have employed w^, /cf = {h — Xi Pi)" and wf^, = ga'~'^i''Jia'- The leading term </)*(g)crp(8)(^ containes both leading- 
twist (LT) and next-to- leading- twist (NLT) contributions in accord with the structures of ap, the terms proportional 



to f^ or jii. The f^ terms would project coUinear qq pair from the meson, while jii terms would not diminish only 
when the qq pair carry noncoUinear momenta. The second step is to substitute the leading parton amplitude '5^ into 
the convolution integral with the meson wave function to extract high twist contributions 

(/)* (g) CTp ® (/i = (/>o ® CTp (8) (/lo + </'o ® ^P ® '/'i + '/'i ® ^P ® <^o + • • ■ , (5) 

where 0o and 0i denote the leading and next-to-leading meson DAs, respectively. The high twist DA (\)\ contains both 
short distance and long distance contributions. The short distance contributions of 0i arise from the noncoUinear 
components of fc^. By the equations of motion, the noncoUinear components of fc^ will induce one quark-gluon vertex 
i7a and one special propagator if^jIxiQ pO|| . Because that the special propagator is not propagating, the quark-gluon 
vertex and the special propagator should be included into the hard function, CTp. In this way, we may factorize 0i as 
(/)i K, (0{^)qwJ^, (0f )" and absorb the short distance piece {<i>x)a into It p. It leads to the third step 

(^l®Gp®(i>X^<^l®ap® ((0f )a • <' ('/'f )"') 

= </.S0(^p.0f)a®<.(</>f)"' , (6) 

where (0f)" containing covariant derivative Z?" —id"' — gA" is implied. Notice that the light-cone gauge w^ • A = 
assures wf^iA" — A". This is legal since each part of Eq. (||) is separately gauge invariant. The second term 
of Eq.(||) also contribute to high twist corrections, as it convolutes with 0o- The momentum factor fc" will be 
absorbed by c/jg to become a coordinate derivative, denoted as k°'(f)Q = (t>i g. Consider another NLT contributions 

fTi w (/fg (g) (CTp)a (Xi w",cj)^ j^ from Fig. 1(b) and (c), where 0"^ containes gauge fields. We have employed the 
approximation that {Wp)^ ® 4>i a '^^ ^^^ leading term of cti. This comes to the fourth step: 

0S ® i^p)c. ® K'(f>i,d + 0S ® (^p)" ^ K4i[a = 't>l® (o^p)i ® ((/>f ) , (7) 

where we have employed (t^i a ^ 't'l A — (0f )"■ However, it will be found that (ap)i diminishes as convoluting with 
twist-4 DA (j)\ (see below definition). Up to NLO in Q~^ , we may drop the {ap)\ term and arrive at the factorization 
for tree amplitudes 

(To + CTl ~ 05 «) O^p «) 00 + 0S «) (ffp • 0f ) «) 01 -f 01 «) ((0t^) • O^p) ® 00 (8) 

where 0i means 0f . There involves only one NLT DA 0i for NLO power corrections. 

To proceed, we need to consider the factorizations of the spin indices, the color indices and the momentum integrals 
over loop partons. For factorization of spin indices, we employ the expansion of the meson DA into its spin components 
as 



Po,i 



E<ir (9) 



where F means Dirac matrix F — 1,7'^, 7'^ 75, cr'^". The factorization of the color indices take the convention that 
the color indices of the parton amplitudes are extract and attributed to the meson DAs. The factorization of the 
momentum integral is performed by making use of the fact that the leading parton amplitudes depend only on the 
momentum fraction variables Xi. The indentity can always be used 

1 
dxi5{xi — ki ■ vl) = I . (10) 

The choice of the lowest twist components 0q ^ of 0o,i is made by employing the power counting. Consider tt meson 
whose high twist DA (^^i'Mi?;ai-aB j^ag ^j^i fcrmion index F and the boson index B. The fcrmion index F arise from 
the spin index factorization for 2F fermion lines connecting DA and parton amplitude and the boson index B denotes 
the no power of momenta in previous coUinear expansion and the no giuon lines as _B = njj + no- We may write 

i 

where A denotes a small scale associated with DA. Spin polarizers e^ denote the combinnation of vectors w^, v,^^ and 
7^. Variable r^ represents the twist of DA 0'. The restrictions over projector g^i'A'-F.ai-aB ^^.^ 



-L/r-"""^-""-""=o (12) 

which are due to the fact that polarizers e^ are always projected by w"^, . The dimension of (^mi ■M-Fiai-aij jg determined 
by dimensional analysis 

d{4>) ^ 3F + B - I (13) 

By equating the dimensions of both sides of Eq. ^V) , one can derive the minimum of r^ 

= 2F + B+l[l-i-lf]. (14) 



T, 



It is obvious from Eq.(n4) that there are only finite numbers of fermion lines, gluon lines and derivatives contributes 
to a given power of 1/Q^- 

B. Collinear Expansion for Arbitrary Loop Orders 

The extension of the collinear expansion for tree diagrams can be straightforward to diagrams containing arbitrary 
loop corrections. The starting point is to notice that the collinear expansion for the one loop corrections in the 
collinear region of the radiative gluons can be written down as 



i=0 j=0 



(15) 



where superscript (0), (1) denote tree and one loop corrections, repsectively. This is because, as the collinear gluons 
with momentum I ~ {Q,\^/Q,X) go through the fermion lines, the valence fermion momenta behave similarly to 
those in the tree level expansion. This leads to the fact that the collinear expansions for one loop amplitudes in 
collinear region can be performed just like for tree amplitudes. The soft gluon corrections can not affect the collinear 
expansion. The cancellations of double logarithms are assured in light-cone gauge by adding ladder and selfenergy 
diagarms. The one loop corrected parton amplitudes are determined by subtracting the amplitudes in collinear and 
soft regions from the full one loop amplitudes . Following the standard considerations |2l[| , the LT paron amplitude 

Wp and NLT paron amplitude ap • {(f>^Y'^' are infrared finite, and the soft divergences are absorbed by 4>q and 
(j>\ . The one loop factorization is derived up to NLT order 

j=o j=o j=o 

j=a i=0 j=0 k=0,l 

+E^o^'^]«EE('/'i'')^^"^'^'^^']«[E ^i'^] ■ (16) 

j=0 i=0 j=0 k=0,l 

The generalization to arbitrary loop orders can be obtained by iteration. Suppose that 

cr = 00 ® (crp)o ® 00 + 0S (^ (o-p)l ® 01 + 0t ® (crp)i ® 00 (17) 

where 



where 



N N N 

•8=0 i=0 i=0 

N 

K)i = Y.^a,t , (18) 

i=0 



TV 






K)i'' = E4^'^-(^f )''"'■'■ (19) 

The above factorization still holds for A^ + 1 order corrections, since the coUinear gluons cannot attach to the parton 
amplitudes. The remaining proof of factorization requires the cancellations of double logarithms of soft divergences, 
the single soft logarithms absorbed by pion DAs and the infradred finitness of the parton amplitudes. This can be 
achieved by standard analysis (see e.g. |^) and it is left to other pubHsh pl[ . 

III. 0(Q"4) CONTRIBUTIONS 

The leading twist contribution of the pion form factor Ft^{Q'^) is expressed as 

(t){x2 )(t>{xi) 

ia;2 
where Q^ = (1 — cos9)Q'^/2 has been employed . Applying (j){xi) — 3f-^Xi{l — a;i)/\/2 into Eq. (pO|), one can get 

tf-"(Q-)= '^^°-^'?'«'^- . (21) 

We use Jtt = 93 MeV in this paper. 

The Feynman diagrams displayed in Fig. 2 contribute to 0{Q^'^). The calculation is straightforward and the result 
reads 

P(t=4).^2^ 256na,iQlg) f JG(xi) + G(.ti)(1 - 2xi)] , ^ 

F^ >{Q)^ ^^ j dxidx2(i>{,X2) X2xl{l-xi) + (a^i^^a:2), (22) 

where G and G mean the twist-4 pion distribution amplitudes |1^] . Note that we have employed the effective coupling 
as{Q^rA with argument Q^g- =< xiX2Q'^ >, where the bracket means the average under the nonperturbative QCD 
vaccum. This is contrary to the conventional treatment that as{Q^) is running with Q^. The scale xiX2Q^ of a^ is 
from the momentum of the exchange gluon. From the factorization point of view, we may write < xiX2Q^ > into 
< xi >< X2 > Q^ t ^y setting a factorization scale Qf^ct to marker the border between the perturbative and the 
nonperturbative dynamics. The scattering angle dependence has been absorbed by Q as Q^ = Q^(l — cos6')^/2(l + 
cos^). The dependence in 9 is irrelevant for power correction, because it is hard to be identified by experiment. 
The contribution from G{xi) is suppressed by the factor (1 — 2xi). Therefore, the contribution mainly comes from 
G(xi) = 3V27i'^/^a;i(l — xi) |lq|. After substituting the distribution amplitudes and performing the integrations, 
there remaines an infrared divergence ||l2| , p2| 

F(-)(g^)^-^^^:^^!iSll:^/^. (23) 



Xl 

It can not be completely resolved under perturbation theory. There requires a resummation over the soft divergences 
as the virtual quark lines become on-shell. One also needs to introduce a jet function to absorb these divergences. 
We will skip the details of the perturbative behavior of these divergences p2|- [l4|j2^ . The initial function for such a 
jet function is of nonperturbative. We denote the effective J function as J{Q2(r) and rewrite the NLO form factor as 



Adding Fi*-^)(Q2) and F^'-'^^Q'^), we obtain 

FAQ') = ^?^^(i Q2 ) + oiQ-') . (25) 

Inspired by the theoretical resuh, we can perform a least x^ (Xmin ~ 7.96742) fit for the data M-M to obtain 



.Fit^^2N _ 0.46895 0.3009^ 



FriQ') = -7^^(l - -7^) • (26) 



The x^ analysis for the data is shown in Fig. 3. It is obvious that the data point at 10 GeV^ is out of allowed errors. By 
igoring the 10 GeV^ data point, we can find that the Q^ — > oo limit of Q^Ft^{Q^) approaches a constant, supporting 
our taking as(Q^rT-) in the above calculations. Comparing the fit formula and the theory formula for Ft^{Q'^), we are 
led to the following conclusions: 

1. The argument Q^ of as{Q^) should be interpreted as an effective Q^rr =< xi >< X2 > Q^ ^- That is we 
need to take asiQ"^) as an effective coupling constant. The average fractions < xi >«< X2 >~ 0.57 and the 
factorization scale Q? , = 1 GeV for Aqqj) = 0.3 GeV. The values oi < Xi > and Q? , depend on the model 
of the pion wave function we have employed (the AS model). If we perform similar analysis by employing CZ 
model 0(x) = 15/^a;(l - a;)(l - 2a;)VV2 and G{x) ^ 15V2tt'^ fi^x{l - x){l - 2xf ||l|l, we may then obtain 
Q2^^^ = 13.12 GeV^ for < x, >w 0.5 or Q\^^^ = 328 GeV^ for < x, >w 0.1. It is clear that CZ model is less 
consistent with PQCD than AS model. 

2. A less model dependent property of the effective coupling constant can be derscribed: The change in Aqq^ 
would affect the location of the average fraction variable < a: > for a fixed factorization scale Qeflt- On the other 
hand, for a fixed < a; >, Qgfj- would vary with Aqqq. Nevertheless, there are only finite possible consistent 
solutions for Aqqq, Qgg- and < a; > can be derived. 

3. The effective value of J function at Q^a is about 0.22. 



4. If we combine the result F~fT^{Q ) of [18 



2/ 8^V.' 

7tV"^ ) — 7^U Q 



F,AQ')^^i^-^4^), (27) 



we may extract an effective running coupling constant uj 

The factor one fourth in a'^ {Q^/4) is from < Xi >w 0.5. The comparison between aj^{Q^) and the usual 
running coupling constant as(Q^) with Aqqq = 0.3 GeV is shown in Fig. 4. One can see that their difference is 
very small for all range Q^ > 1 GeV^. To reduce their differences, higher order power corrections are required. 
The effective running coupling constant af /tt is smaller than 0.2 for Q^ > 1 GeV^. At the accuarcy we work, 
it directly showes that PQCD is applicable to exclusive processes for momentum transfer of Q^ > 1 GeV^ 
{a,/TT < 0.5). 

IV. DISCUSSIONS 

From the above analysis, we may draw the conclusion that PQCD is applicable to 7*7r — > tt at the energy region 
accessable by the experiment. In our analysis, we only employed the simplest models without invoking any nonpertur- 
bative arguements. From our result, we may try to understand the evidence showed by Isgure and Llewenlly Smith. 
In their analysis, they employed the fraction function 

fie)=f\xf\yOixy-e)^t(-}m (29) 

Jo Jo ^y 



to evaluate the percentage of perturbation contributions. The parameter e describes a cut-ofF on xy in order to keep 
higher-twist contributions and higher-order effects at a reasonable small level. Since /(e) is a mixing of perturbative 
and nonperturbative contributions, the perturbation parts are those corresponding to large values of xyQ'^ . For the 
AS model, /(e) reaches 90 % for e = 1/150. It implies that the naive perturbative contribution is 90% legal for 
Q^ — 150 GeV^. It is noted that they employed a common setting for the scale of the running coupling constant 
cts{Q^) in the pion factor factor, the small value of as{Q'^) for Q^ > 1 GcV^ leads to their conclusion that PQCD is 
inapplicable to exclusive process. 

However, the real value of the coupling constant or the factorization scale involving in process 7*7r — > tt can not be 
determined a priori by perturbation theory. On the contrary, it should be determined from other methods, e.g. the 
experiment. In fact, /(e) is reasonable only for the asymptotic region at Q^ -^ oo. As shown above, it is indeed. But 
the lack of experiment in such a high Q^ region, it is difficult to determine the factorization scale. The analysis of 
Isgur and Llewellyn Smith employed the leading result, which should be valid for Q^ — > cxd. Therefore, they needed to 
extrapolate the leading result from high Q^ to low Q^. It seems that the conclusion of Isgur and Llewellyn Smith is 
weak. Because we have employed the leading and sub-leading contributions in the pion form factor, we can control the 
Q^ behaviors of pion form factor for both high and low Q^. Therefore, our conclusion is stronger and more reliable. 

The power corrections to the pion form factor have also been calculated by employing the projection for the pion 

P{k,p-k) = i75{l + ^^^^P.'^'-P'. + ^^^^P^'y^-k^.} (30) 

where k is the momentum carried by the valence quark of the pion, x represents the fraction x — k ■ p' /p ■ p' and p 
and p' denote the momenta carried by the initial and final state pions. The result appears as 



Q2 ^^ ' Q2^2 



FAQ') = ^^f^^{i + 7^2^:2 J'iQ')} ■ (31) 



The function J{Q^) represents the J function introduced before. The factor ^^ > 1.4 GeV results in a violation of 
perturbation principle for the region of Q^ being 2 ~ 15 GeV^. The positive sign in 0((5^^) power contributions 
seems difhcult to expain the x^ fit formula of the data. 

Another factor that would affect the leading result may arise from the contributions of Sudakov form factors. 
However, as shown in Ul that the Sudakov form factor falls fastly than any power suppression for low Q^ . Then, the 
low Q^ behaviors of the pion form factor is mainly controlled by the power corrections (the sub-leading contributions). 

V. CONCLUSIONS 

We have developed a power expansion scheme for 7*7r —^ tt. This expansion scheme also has been shown to be 
compatible with QCD factorization. The NLO power corrections to the pion form factor have been calculated. With 
the help of the formula of the pion form factor upto Q^**, we may perform a least x^ analysis for the data. By 
comparing the x^ fitting formula and the theory formula of the pion form factor, we arrive at the conclusion that, 
due to the nonperturbative QCD vaccum, the strong coupling constant in the theory formula should be identified as 
an effective coupling with factorization scale equal to 1 GeV. In addition, the averaged fraction variable locates at 0.5 
in consistency with the AS model for pion wave function. The CZ model is fail to give a consistent explanantion for 
the data, because its relative factorization scale is in the range of 3 ^ 18 GcV. 

We also derive the Q^ behavior of the effective coupling, the effective running coupling constant. From the Q^ 
analysis of the effective running coupling constant, we may prove that PQCD is applicable to exclusive process for 
momentum transfer larger than 1 GeV^. This is a direct proof showing that PQCD is applicable to exclusive process 
for Q^ > 1 GeV^. In addtion, the analysis performed by Isgur and Llewellyn Smith showing that PQCD is inapplicable 
to exclusive process, can also be undertood that the leading contribution is valid in the asymptotic region and can 
not be applied for low Q^. 

The result obtained in this paper can be easily extended to other hard processes involving two light mesons, such 
a.s r], p ^ 7r7, etc. 

From the result of this paper and pq ], we have a suggestion that the power corrections in exclusive processes 
are sizable and require detail investigations. Over past decades, we have accumulated abundant data for exclusive 
processes. Most data are in low energy region. This is the place where the power correction become important and 
can not be negligible. The leading or the asymptotic contribution only gives information about the hadron wave 



function, the static property of QCD, while the power correction can reveal the dynamics. This is why the power 
correction plays an important role in analysis of data. 
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Figure Caption 

Fig.l The leading order diagrams for j*tt —^ tt. 

Fig. 2 The next-to-leading-twist (NLT) diagrams for j*tt -^ tt. The propagator with one bar means the special 
propagator. 

Fig.3 Plot of the_ least x^ fit (solid hne) and C.L.= 99.73% (dash hne) for Q^F^{Q'^). The experimental data 
are taken from 



Fig. 4 The comparison between the effective running coupling constant (solid line) and the perturbative running 
coupling constant (dash line). 
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